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Abstract 

We report photometric and radial velocity observations of the XO-4 transiting 
planetary system, conducted with the FLWO 1.2m telescope and the 8.2m Subaru 
Telescope. Based on the new light curves, the refined transit ephemeris of XO-4b is 
P = 4.1250828±0.0000040 days and T c [BJD TDB ] =2454485.93323±0.00039. We mea- 
sured the Rossiter-McLaughlin effect of XO-4b and estimated the sky-projected angle 
between the stellar spin axis and the planetary orbital axis tobeA = -46.7°li;lo. This 
measurement of A is less robust than in some other cases because the impact param- 
eter of the transit is small, causing a strong degeneracy between A and the projected 
stellar rotational velocity. Nevertheless, our finding of a spin-orbit misalignment sug- 
gests that the migration process for XO-4b involved few-body dynamics rather than 
interaction with a gaseous disk. In addition, our result conforms with the pattern re- 
ported by Winn et al. (2010, ApJL, 718, L145) that high obliquities are preferentially 
found for stars with effective temperatures hotter than 6250 K. 

Key words: stars: planetary systems: individual (XO-4) — stars: rotation - 
techniques: radial velocities — techniques: spectroscopic — techniques: photometric 



Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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1. Introduction 



Observations of the Rossiter-McLaughlin effect (hereafter the RM effect: Rossiter f 924; 
McLaughlin f924; Ohta et al. 2005; Gaudi & Winn 2007; Hirano et al. 2010) for transiting 
planetary systems enable us to measure the sky-projected angle between the stellar spin axis and 
the planetary orbital axis (also called the projected spin-orbit angle, or the projected obliquity). 
Previous measurements of this angle have revealed the existence of well-aligned prograde orbits 
(|A| < 10°; see, e.g., HD 209458b: Queloz et al. 2000; Winn et al. 2005, HD 189733b: Winn et 
al. 2006), highly-inclined prograde orbits (10° < |A| < 90°; XO-3b: Hebrard et al. 2008; Winn 
et al. 2009b, HD80606b: Moutou et al. 2009; Winn et al. 2009c; Hebrard et al. 2010, WASP- 
14b: Johnson et al. 2009, CoRoT-lb: Pont et al. 2010), and retrograde orbits (90° < |A| < 180°; 
HAT-P-7b: Narita et al. 2009a; Winn et al. 2009a, WASP-17b: Anderson et al. 2010; Triaud et 
al. 2010, WASP-33b: Cameron et al. 2010, WASP-8b: Queloz et al. 2010, WASP-15b: Triaud 
et al. 2010). 

The discovery of highly tilted orbits, and retrograde orbits, has been taken as evidence 
in favor of "few-body" theories for planetary migration, such as planet-planet scattering and 
subsequent tidal evolution (planet-planet scattering models, e.g., Rasio & Ford 1996; Marzari 
& Weidenschilling 2002; Nagasawa et al. 2008; Chatterjee et al. 2008), or the Kozai mechanism 
(Kozai 1962) caused by a distant companion and subsequent tidal evolution (Kozai migration 
models, e.g., Wu & Murray 2003; Takeda & Rasio 2005; Fabrycky & Tremaine 2007; Wu et al. 
2007). In this way, the RM effect has already played an important role for studies of planetary 
migration processes. Further measurements of the RM effect are important in order to constrain 
theoretical models by learning the statistical distribution of spin-orbit alignment angles, and 
by discriminating planetary migration mechanisms for individual planetary systems. For this 
purpose, we have continued measurements of the RM effect and have also initiated direct 
imaging observations, with the Subaru 8.2m telescope, of systems with misaligned planets to 
search for evidence of perturbing stellar companions (e.g., Narita et al. 2007; Narita et al. 
2010b). 

The subject of this letter is the transiting exoplanet XO-4b, discovered by McCullough 
et al. (2008) (hereafter MC08). The host star XO-4 is a relatively bright (V = 10.7) F5 star, 
and XO-4b is a giant planet (M p = 1.72 Mj up , R p = 1.34 -Rj up ; MC08) with an orbital period 
of P = 4.125 days. We used the FLWO 1.2m telescope for photometric observations and the 
Subaru 8.2m telescope for radial velocity (RV) observations. Due to the rapid rotation of XO-4 
(Vsinl s = 8.8 ±0.5 km s _1 ; MC08), the amplitude of the RM effect for XO-4b was expected to 
be large. In addition, recently Winn et al. (2010) reported that hot stars (T e g- > 6250 K) with 
hot Jupiters tend to have high obliquities. Since XO-4 is a hot star (T e g = 6397 ± 70 K; MC08) 
with a hot Jupiter, we may use this system as a further test of the proposed pattern. We 
report our photometric and RV observations in section 2, and describe our analysis procedures 
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in section 3. We present our results and discussions in section 4. Finally, we summarize the 
findings of this letter in section 5. 

2. Observations 

2.1. Light Curves 

Photometric observations were conducted with Keplercam on the 1.2m telescope at the 
Fred Lawrence Whipple Observatory (FLWO) on Mount Hopkins, Arizona. Three partial 
transits, and one complete transit, were observed with a Sloan Digital Sky Survey (SDSS) z- 
band filter. The observing dates were UT 2008 October 5 (egress), 2008 October 13 (ingress), 
2009 September 8 (ingress), and 2010 February 8 (complete). KeplerCam is equipped with a 
4096x4096 CCD, giving a square field of view 23.1' on a side. We used 2x2 binning, giving a 
scale of 0.68" per binned pixel, a readout and setup time of 11 s, and a typical readout noise of 7 
e~ per binned pixel. The cadence was typically 2 min for the observation on UT 2010 February 
8, and about 50 s for the other observations. Standard IRAF 1 procedures were applied to the 
observed frames for bias subtraction, flat-field division, and aperture photometry. The flux of 
XO-4 was divided by the sum of the fluxes of 11 nearby reference stars in the field of view. To 
account for systematic effects due to differential airmass extinction and imperfect flat-fielding, 
we corrected the baseline of the observed flux as a function of time and airmass. Photometric 
uncertainties were temporarily set equal to the standard deviation of the out-of-transit relative 
flux on that night. As explained in subsection 3.1, we subsequently rescaled these uncertainties 
to account for time-correlated noise. The upper panel of figure 1 shows the composite transit 
light curves. 

2.2. Radial Velocities 

We observed XO-4 on UT 2010 January 14, 2010 January 15, 2010 February 4, and 2010 
May 17 with the High Dispersion Spectrograph (HDS: Noguchi et al. 2002) on the 8.2m Subaru 
Telescope at Mauna Kea, Hawaii. The observations on 2010 January 15 covered a full transit of 
XO-4b. We employed the standard I2a setup of the HDS, covering the wavelength range 4940 A 
< A < 6180 A and a slit width of 0"4, corresponding to a spectral resolution of about 90000. 
We used the iodine gas absorption cell for precise differential RV measurements. The exposure 
times were 15-20 min, yielding a typical signal-to-noise ratio (SNR) of approximately 80-100 
per pixel. We processed the observed frames with standard IRAF procedures and extracted 
one- dimensional spectra. We computed relative RVs following the algorithm of Butler et al. 
(1996) and Sato et al. (2002). We estimated the internal error of each RV based on the scatter 
among the RV solutions from ~4 A segments of each spectrum. The typical internal errors were 

1 The Image Reduction and Analysis Facility (IRAF) is distributed by the U.S. National Optical Astronomy 
Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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9-15 m s _1 . As was the case for TrES-4 (Narita et al. 2010a), the internal errors are larger than 
some other cases of spectra with a similar SNR due to the star's broad absorption features. 
The observed RVs and internal errors are given in table 1. We note that we adopt BJDtdb, 
the Barycentric Julian Date in the Barycentric Dynamical Time standard, for all time stamps 
in this letter, as advocated by Eastman et al. (2010). 

3. Analyses 

3.1. Refined Transit Ephemeris 

First we fitted the FLWO transit light curves in order to determine appropriate data 
weights for our subsequent analysis, and to refine the transit ephemeris of XO-4b. We used 
the analytic formula for transit light curves given by Ohta et al. (2009). The free parameters 
were the midtransit time T c for each transit, the ratio of radii of the planet and star R p /R s , 
the orbital inclination i, and the orbital distance in units of the stellar radius a/R s . We 
here adopted the orbital period P = 4.12502 days, the origin of mid-transit time of XO-4b 
T c (0) = 2454485.93295 ±0.00040 in BJD TD b, the stellar mass M s = 1.32 M reported by MC08. 
We held fixed one of quadratic limb-darkening coefficients U\ at the value 0.13, based on tables 
of Claret (2004), and allowed 112 to be a free parameter. 

We determined the optimal parameter values by minimizing the x 2 statistic, 



using the AMOEBA algorithm (Press et al. 1992). In this equation, /j )0 b s are the relative flux 
data points, and /j jCa ic are the values calculated based on the analytic formulae and a particular 
choice of model parameters. 

We then computed the residuals for each dataset, and averaged the residuals into M 
bins of N points. We calculated the standard deviation of the binned data, 0"N, o bs and 



where o\ is the standard deviation of the residuals for each dataset. To account for increased 
uncertainties due to time-correlated noise (so-called red noise; see e.g., Gillon et al. 2006; Winn 
et al. 2008), we computed a red noise factor j3 = crN,obs/o"N,idcai for various N (corresponding 
to 10-20 min), and multiplied the photometric uncertainties of each dataset by the maximum 
value of p. We estimated la uncertainties of free parameters based on the criterion A% 2 = 1.0. 
Table 2 summarizes the rms, the red noise factor, and the mid-transit time of each transit light 
curve. One important result is that i = 90° is allowed within the la level. This was problematic 
for subsequent RM modeling, because the two key parameters describing the RM effect (A and 
Vsinl s , see below) are strongly degenerate when i = 90°. 

We then fitted all the observed mid-transit times with a linear function 



x 2 = E 




(i) 




(2) 



4 



T c (E) = EP + T c (0), (3) 

where E is an integer, in order to refine the orbital period of X0-4b. We note that we added 
0.00075 to the value of T c (0) reported by MC08 to convert the time standard from HJDutc 
(the Heliocentric Julian Date in the Coordinated Universal Time) to BJDtdb (Eastman et al. 
2010). 

The refined transit ephemeris for XO-4b is P = 4.1250828 ± 0.0000040 days and T c (0) = 
2454485.93323 ±0.00039 in the BJDtdb system. We note that the linear fit gave \ 2 = 6.99 and 
all O-C residuals were within 1.5<t uncertainties of T C (E). The uncertainty in the predicted 
time of the transit of UT 2010 January 15 (the night of our RM measurement) is about 1 min, 
which is small enough to be of no concern for our RM modeling. 

3.2. RM Model and Joint Fitting 

To describe the RM effect of XO-4b, we derived a formula based on the procedure 
described in Winn et al. (2005), Narita et al. (2009b), and Hirano et al. (2010). We assumed a 
rotational broadening kernel of Vsinl s = 8.8 km s _1 (MC08) and the quadratic limb-darkening 
parameters u\ = 0.35 and U2 = 0.35 for the iodine absorption band based on the tables of Claret 
(2004). The derived formula was 



Av 



-fv P 



1.6159-0.83778 



(4) 



s V sin I 

In this formula, Av is the anomalous radial velocity due to the RM effect, / is the loss of light, 
and v p is the subplanet velocity (the radial component of the stellar rotation velocity at the 
position being hidden by the planet). 

We simultaneously fitted the Subaru RVs and the FLWO light curves. We fixed P and 
T c (0) at the central values reported above. We did not fit the RVs reported by MC08 because 
the errors of the MC08 RVs are much larger than those of the Subaru RVs. Thus, those 
data have little impact on the results. In addition to the free parameters for the photometric 
data, we added free parameters for the radial velocity semiamplitude K, the sky-projected 
angle between the stellar spin axis and the planetary orbital axis A, the sky-projected stellar 
rotational velocity Vsml s , and the overall velocity offset of the Subaru dataset v±. Note that 
we fixed the eccentricity e to zero, and the argument of periastron u was not considered. (When 
e and uj were included as free parameters, the eccentricity was found to be consistent with zero 
within la.) The \ 2 statistic for the joint fit was 



X 



E 



f j.obs 



i,calc 



E 



^j',obs ^j',calc 



(J, 



+ 



Vsinis-8.8 
05 



where t^obs were the observed RVs and Vj >ca x c were RVs calculated based on the Keplerian 
motion and on the RM formula given above. The last term was an a priori constraint on 
U sin J, to require consistency with the value reported by MC08 within its uncertainty. This 
term was necessary to obtain a meaningful conclusion about A, as a consequence of the low 
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impact parameter (near-90° inclination) of the transit (see subsection 3.1.). 
4. Results and Discussion 

The lower panel of figure 1 shows the RVs around the transit phase, along with the 
best-fitting model (the solid line). For reference, we have also plotted a model RV curve under 
the assumption of perfect spin-orbit alignment (A = 0°, the dotted line). We found small 
asymmetry of the RM anomaly with respect to the midtransit time, suggesting that the transit 
is not exactly equatorial and that there is a significant spin-orbit misalignment. 

Figure 2 plots the full RV dataset in BJDtdb (left panel) and in phase (right panel). No 
significant long term trend was observed in the residuals. Table 3 summarizes the optimal values 
and uncertainties for the system parameters. As shown in table 3, the free parameters other 
than the RM-related parameters (A and Vsin/ S ) are consistent with those reported previously. 
For the RM-related parameters, we found A = — 46.7°lg;{o and Vsinl s = 8.9 ±0.5 km s^ 1 . We 
also computed the 3a confidence level of A by A% 2 = 9.0 and found A = — 46.7° li^g - Figure 3 
shows a x 2 contour map in X-Vsml s space. The solid lines represent contours for A% 2 = 1.0, 
Ax 2 = 4.0, and A% 2 = 9.0 from the inside to the outside. The orbital axis of XO-4b appears to 
be tilted relative to the rotational axis of XO-4 in sky projection at about the 3a level. 

It should be noted that the confidence with which A is found to be nonzero depends 
critically on the a priori constraint that we applied on V sin/,.. This is because when the transit 
impact parameter is small (or equivalently when the orbital inclination is near 90°), there is a 
strong degeneracy between ^sin/ s and A. In such cases, the RV data alone are only able to 
provide good constraints on the quantity ^sin/ s cosA, while ^sin/ s sinA is poorly constrained. 
When we tried the joint fitting without the a priori constraint, we found \^sin/ s cosA = 8.4 ± 
1.1 km s" 1 , Vsml s = 10.0^ 2 2 ° km s _1 , A = -54.6° and z = 89.1° ±0.9° at the la level. 
In this case the confidence in the spin-orbit misalignment was lowered, and values of Vsinl s 
differing strongly from that reported by MC08 were allowed within la. 

As long as the spectroscopic determination of Vsinls is valid, our result of a strong 
spin-orbit misalignment is secure. However, in order to make the result more robust, and less 
dependent on the spectroscopic determination of Vsinl s , it would be necessary to obtain a 
more stringent constraint on the inclination i. The most straightforward way to do so would be 
to obtain a more precise light curve. It is difficult to observe complete transits from the ground, 
because of the comparatively long orbital period and long duration of the transits. However, 
the system is in the continuous viewing zone of the Hubble Space Telescope (HST), as noted 
by MC08. A complete transit light curve based on HST observations is therefore the highest 
priority for pinning down the spin-orbit angle. 

From a theoretical point of view, it is interesting that the planetary orbit is highly 
tilted, given that the orbital eccentricity of XO-4b is consistent with zero at this point in time. 
However, the current bounds on the orbital eccentricity are coarse; improving the determination 
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of this parameter is also a priority for future work. 

Recently, Winn et al. (2010) pointed out that hot Jupiters around relatively hot (T c g- > 
6250 K) stars preferentially have tilted orbits. According to the empirical rule, XO-4, which 
has an effective temperature of T e g = 6397 ±70 K (MC08), would have a tilted orbit. Our 
results agree with this prediction. 

In addition, Schlaufman (2010) has recently sought evidence for spin-orbit misalignments 
along the line of sight, based on the comparison of the measured V sin Is and the plausible 
range of rotation velocities for a star of the given mass, radius, and age. He argued that if a 
value = (V sml s sim — V sin I s>0 bs) / J ^sim + a obs * s l ar g er than 2.9 (a threshold based on the 
SPOCS catalog: Valenti & Fischer 2005), the stellar rotational axis is very likely to be inclined 
with respect to the line of sight. For a transiting planet host, this would imply a spin-orbit 
misalignment. (He also found that hot stars were more likely to have high obliquities, although 
he phrased the correlation in terms of stellar mass rather than effective temperature.) For 
XO-4, the value is O = 2.61, which is only just below the (somewhat arbitrary) threshold of 
2.9. This suggests that XO-4 has a smaller projected rotational velocity than expected, and 
can be taken as supporting evidence for a spin-orbit misalignment in the XO-4 system. 

5. Summary 

We conducted photometric transit observations with the FLWO 1.2m telescope, and 
RV observations (including an RM measurement) with the 8.2m Subaru Telescope, for the 
exoplanet XO-4b. We refined the transit ephemeris for XO-4b based on the new transit light 
curves, and modeled the RM effect of XO-4b by jointly fitting the photometry and RVs. In 
conjunction with the previous spectroscopic estimate of the projected stellar rotational velocity 
Vsinl s , we found evidence for spin-orbit misalignment in this system. In this case the RM 
data do not provide an strong independent constraint on Vsml s , because of the low impact 
parameter of the transit. To complement the RM measurement in this system, a more stringent 
constraint on the orbital inclination would be useful to break the fitting degeneracy between A 
and \^sin/ s . 

The spin-orbit misalignment of XO-4b suggests that this planet migrated through a 
mechanism that excited its orbital inclination, although the eccentricity of XO-4b seemed to be 
quite low. Furthermore, the spin-orbit misalignment of XO-4b is consistent with the findings by 
Winn et al. (2010) and Schlaufman (2010) that hot stars with hot Jupiters have high obliquities. 
Further follow-up observations to determine a more accurate spin-orbit alignment angle in this 
system would be still interesting. 

This letter is based on data collected at Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. We acknowledge the support for our Subaru 
HDS observations by Akito Tajitsu, a support scientist for the Subaru HDS. The data analysis 
was in part carried out on common use data analysis computer system at the Astronomy Data 
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Fig. 1. Top panel: Phased transit light curves taken with the FLWO 1.2m telescope. Four different 
epoch data are combined. The solid line represents the best-fit curve based on the optimal free parameters 
listed in table 3. Second panel: Residuals of the photometric data from the best-fit model. Third panel: 
Subaru radial velocities and the best-fit curve (the solid line) as well as a model curve which assumed 
A = 0° (the dotted line). Bottom panel: Residuals of radial velocities from the best- fit model. 




Fig. 2. Upper panels: Subaru radial velocities in BJDtdb (left panel) and in phase (right panel). Lower 
panels: Residuals from the best-fit model. No significant long term trend is apparent. 
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Fig. 3. A x 2 contour map in X-V sin/ s space. The solid lines represent contours for A\ 2 = 1.0, Ax 2 = 4.0, 
and Ax 2 = 9.0 from the inside to the outside. 
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Table 1. Radial velocities obtained with the Subaru/HDS. 
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2455212.04223 


-35.98 


9.33 


2455212.05343 


-29.42 


9.41 


2455212.07586 


-18.34 


11.22 


2455212.08708 


-42.78 


9.52 


A H r" f AAA f~\ /~* A O /~* 

2455232.06436 


114.02 


18.42 


2455232.07383 


97.09 


13.09 


2455232.08330 


107.35 


12.66 


2455333.73171 


-15.49 


15.33 


2455333.74176 


47.77 


12.39 


2455333.75088 


23.27 


11.23 


2455333.76001 


26.43 


13.70 


2455333.76915 


a£67 


13.64 



Table 2. Summary of the XO-4 transit light curves. 



UT Date 


iLpocn 


KMb 


Red noise factor 


Mid-transit time 


Note 






h. 


1U a 


o 

p 


i3jL»TDB 




2008 Jan. 


20 









2454485.93295 ± 0.00040 


MC08 


2008 Oct. 


05 


63 


12.7 


1.293 


2454745.81494 ±0.00105 


Egress 


2008 Oct. 


13 


G4 


15.4 


1.497 


2454749.93969 ±0.00114 


Ingress 


2009 Sep. 


08 


144 


17.2 


1.002 


2455079.94672 ±0.00110 


Ingress 


2010 Feb. 


08 


182 


8.2 


1.807 


2455236.69753 ±0.00065 


Complete 
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Table 3. System parameters of XO-4 derived from the joint fitting. 



Parameter 



Value 



Uncertainty 



Note 



Adopted parameters 



P [days] 

T c (0) [BJD TDB ] 

M S [M ] 

Vsinl s [km s _1 ] 

m 

e 



4.1250828 ±0.0000040 this work 

2454485.93323 ±0.00039 this work 

1.32 ± 0.02 MC08 

8.8 ±0.5 MC08 

0.13 Claret (2004) 

assumed 



Best-fit parameters with the a priori constraint 

± 6.2 
±0.5 



K [m s" 1 ] 
V svo-Ig [km s -1 ] 

Ml 

a/R s 
Rp/R s 

V\ [m s _1 ] 
rms [m s _1 ] 
X 2 /v (HV) 
X 2 lv (LC) 



168.6 
8.9 
-46.7 
7.68 
0.0881 
88.8 
0.35 
-0.1 
11.50 
31.63/37 
487.22/906 



+8.1 
-6.1 

± 0.11 

± 0.0007 
± 0.6 
± 0.11 
± 2.9 



±S3 (3-) 



Derived planet parameters 
M p [M Jup ] 1.78 
Rp[Rju P ] 1-33 



± 0.08 
± 0.05 



this work 
this work 



Comparison with previous literature 

K [m s- 1 ] 163 ± 16 MC08 

a/R s 7.7 ± 0.2 MC08 

R p /R s 0.089 ± 0.001 MC08 

i [°] 88.7 ± 1.1 MC08 

M p [M Ju p] 1.72 ± 0.20 MC08 

Rp[Rj up ] 1.34 ± 0.048 MC08 

u 2 0.36 Claret (2004) 
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